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Abstract. This paper presents a new spectral evolutionary model of galaxies, properly taking the effects of nebular 
emission and pre- main sequence evolution into account. The impact of these features in different photometric filters 
is evaluated, along with the influence that variations in the physical conditions of the gas may have on broadband 
colours, line ratios and equivalent widths. Inclusion of nebular emission is demonstrated to radically change the 
predicted ultraviolet, optical and near-infrared colours during active star formation. Pre-main sequence evolution 
is also seen to give a non-negligible contribution to the luminosity in the near-infrared during the first few millions 
years of evolution and should not be omitted when very young systems are being modelled. Finally, we present a 
comparison of our predictions to observations and two other recent codes of evolutionary synthesis. 
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1. Introduction 



Ever since the pioneering work of Tinsley ( |l972| ) , spectral 
evolutionary models (SEMs) based on stellar population 
synthesis have proven to be important tools in the analysis 
of the st ar formation history of galaxies (e.g. Leitherer et 
al. 1996 ). The Sloan survey (e.g. Gunn et al. 1995 ) and 
the Anglo- Australian 2dF survey (e.g. Colless 199?) will 



provide a vast number of galaxy spectra that, with the 
aid of SEMs, may be used for important statistical studies 
of the chemical evolution, initial mass functions and star 
formation histories of galaxies. 

Taking the effects of nebular emission into account be- 
comes essential when deriving the broadband colours of 
galaxies subject to substantial star formation rates. This 
does not only apply to local starbursts, but could also have 
implications for various Hubble types during the first few 
hundred million years of evolution, i.e. the maximum age 
of the most distant galaxies now observed. Despite this 
fact, results are continuously being published where this 
problem is completely overlooked, mainly due to the lack 
of suitable models. It is therefore a necessary and impor- 
tant step to include a proper treatment of nebular emission 
at a wide range of metallicities. 

In connection with our work with metal poor galax- 
ies we have previously developed SEMs from different ap- 
proaches, all of them relevant at low metallicities and 



inclu ding a nebular component (Bergvall 1985 , Olofsson 
1989 ). Here we present a new, more flexible and up-to- 



date model intended for the interpretation of low to in- 
termediate mass galaxies within the full range of ages and 
metallicities observed. Improvements compared to similar 
codes include pre-main sequence evolution and a sophisti- 
cated treatment of nebular emission. 

Since our main interest lies in the evolution of low 
to intermediate mass galaxies, chemical evolution is not 
treated at this stage. Even though a continuously increas- 
ing metallicity is definitely expected for galaxies under a 
closed-box scenario, most galaxies do, most certainly, un- 
dergo infall of gas. Besides hydrogen and helium, various 
amounts of heavy elements should also participate in the 
infall process. Moreover, selective loss of chemical elements 
from stars within a galaxy is also expected, in particu- 
lar in low-mass system (de Young et al. 1990, D'Ercole 
et al. 1999). This effect could even be important in mas- 



Send offprint requests to: 
Erik . Zackrisson@astro .uu.se 



Erik Zackrisson, e-mail 



sive galaxies if the star forming region resides far from 
the centre of the gravitational potential. Under specific 
conditions, the amount of heavy elements may even de- 
crease during certain episodes in the lifetime of a galaxy 
(Olofsson 199!:). A proper treatment of chemical evolu- 
tion for low- to intermediate mass objects would require 
complete knowledge of the dynamics of the system, the 
time-scales of infall, the gas consumption rate and the 
ejection of chemical elements. In lack of a better prescrip- 
tion, we therefore assume constant metallicity evolution 
of each separate stellar population within a galaxy. The 
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spectral energy distributions (SEDs) produced this way 
may then be mixed to simulate the effect of a metallic- 
ity distribution. We claim that this approach should be 
no less appropriate than closed box models subject to ad 
hoc modifications of the rate of infall and mass-loss nec- 
essary to match the metallicities observed. In the future 
we will however investigate the possibility to derive in- 
formation about the chemical evolution directly from a 
tree-SPH model currently being developed at our insti- 
tute (Pharasyn 2001), including a realistic treatment of 
star formation processes. 

The treatment of dust is another important but sim- 
ilarly difficult issue since the size, chemical composition 
and distribution of dust particles may vary considerably 
and lead to large differences in the global spectral distri- 
bution. Simple geometries of dust-star distributions will 
however be included in a future version of the model. 

In this paper, we present the code at its present stage of 
development and investigate the impact of nebular emis- 
sion and pre-main sequence evolution on the modelling of 
the SEDs of galaxies. In section 2, we outline the machin- 
ery on which the code is based - in terms of the compu- 
tational method, stellar evolutionary tracks, stellar atmo- 
spheres, photoionization model and general assumptions 
used. In sections 3 and 4 we investigate the magnitude 
and duration of nebular and pre-main sequence effects in 
different photometric filters and for different metallicities. 
Section 5 includes a comparison of the colours predicted 
by our code to those of two other SEMs. The colours, 
line ratios and equivalent widths predicted are then tested 
against observations in section 6. Section 7 contains a sum- 
mary of our findings and a few concluding remarks. 

2. The model 

The spectral evolutionary synthesis code presented in this 
paper uses a set of stellar evolutionary tracks, theoretical 
stellar atmosphere spectra and a photoionization model to 
predict colours and SEDs for galaxies at ages between 0-15 
Gyr. The stellar component is represented by a single stel- 
lar population evolving at constant metallicity 2=0.001, 
0.004, 0.008, 0.020 or 0.040, whose initial mass function 
(IMF) follows a single-valued power law 

$(m, a) = - — = Cm 
dm 

within the mass interval [m\, m u ] . The normalization con- 
stant C is here determined by the total mass converted 
into stars. The model allows the use of any upper and 
lower mass limit between 0.08-120 M Q and any value of 
a. 

The star formation history can be modelled using a 
constant star formation rate (SFR) of arbitrary duration 
(r): 

Mo 



or, alternatively, a SFR subject to exponential decrease 
over time: 

SFR(*,t) = — exp(-~) 

r r 

where r is the e-folding decay rate of the SFR. In both 
cases Mq represents the total mass available for star for- 
mation. 



2.1. Stellar evolutionary tracks 

In order to follow the evolution of stars between 0.08-120 
M Q at metallicities 2=0.001, 0.004, 0.008, 0.020 and 0.040 
through all important stages of evolution, stellar tracks 
from several different authors have been adopted. 

The tracks developed by the Genev a group 

Schaerer 



1992 



(Charb onnel ct al. |1993| Schallcr ct al. 
et al. 19934 [3 Charbonnel et al. 1996|) are used to 
follow the evolution of stars from the zero age main 
sequence (ZAMS) to the end of the C-burning phase for 
7 < M/Mq < 120 at all Z, to the end of the early asymp- 
totic giant branch (E-AGB) for 2 < M/M Q < 5 and up 
to the He-flash for 0.8 < M/Mq < 2 at 2=0.001, 0.004 
and 0.008. At 2=0.020 and 2=0.040, these tracks follow 
the evolution to the E-AGB throughout the entire mass 
interval 0.8 < M/Mq < 5. Evolution from the He-flash 
to the E-AGB for stars between 0.8 < M/Mq < 2 at 
2=0.001, 0.004, 0.008 is instead covered using the tracks 
by Castellani et al. (1991), since the horizontal branch 
(HB) morphologies generated by the Geneva tracks are 
known to become very red regardless of metallicity. This 
strategy is further discussed in section 2.3. 

Evolution along the pre-main sequence (PMS) has 
then been added for stars b etwee n 0.8 < M/Mq < 7 us- 
ing the tracks by Bernasconi 1996 , and from the E-AGB to 
the white dwarf and planetary nebula stage for stars in the 
range 0.8 < M/Mq < 5 using the hyd rogen - and helium 
burning tracks by Vassiliadis & Wood ( 1994 ). In order to 
make the zero age post-AGB mass and main fuel a unique 
function of ZAMS mass, all tracks computed using artifi- 
cially enhanced or reduced mass-loss rates on the red giant 
branch (RGB) were excluded from the set. This strategy 
will not affect the final predictions since the tracks actu- 
ally employed cover the same range in log(g) and T c g and 
will be interpolated during the course of calculation. 

The gap between E-AGB and the post-AGB has been 
bridged using characteristic luminosities and time-scales 
(~ 10 5 yr) of the thermally pulsing AGB phase (Willson 
1998, 2000) for all stars less massive than 5 M©. These 



SFR(i, t) = 
SFR(t, 



< t < T 



t> T 



numbers are based on a combination of theoretical mass 
loss rates and observations, in which selection effects have 
carefully been taken into account. 

Low-mass tracks from the Lyon-group (Chabrier & 
Baraffe 1997) have been adopted to encompass the evolu- 
tion of stars with masses in the range 0.08 < M/Mq < 0.8. 
Since these tracks terminate at an age of 10 Gyr, it has 
been necessary to make linear extrapolations in log(i) in 
order to obtain a final data point at an age of 15 Gyr. 
These models use non-grey atmospheres, which resolves 
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the temperature-luminosity discrepancies at the bottom 
of the low-mass sequence present in previous generations 
of low-mass models. Kroupa & Tout (1997) compared sev- 
eral low-mass models and concluded that the ones created 
by the Lyon-group provided the best agreement with ob- 
servations. 

All tracks not available at the metallicities used in the 
Geneva set have been been linearly inter- and extrapo- 
lated in log(Z) to cover all five metallicities chosen for 
our model, and, in the case of the post-AGB tracks, also 
to cover the full mass range at Z=0 .001. The ZAMS to 
E-AGB tracks by Charbonnel et al. 



1996 have been ex- 



trapolated up to Z=0.040 using a procedure that closely 
matches each track to the original ZAMS to He-flash 



tracks by Schaerer et al. (1993a) during common stages 
of evolution. 



2.2. Important approximation connected to PMS 
evolution 

PMS stars are assumed not to contribute to the total lu- 
minosity of the galaxy before emerging from their dusty 
parental cocoons. The time-scales of these processes have 



been adopted from Bernasconi (1996). Even though this 
correction is probably appropriate in the optical region, as 
illustrated by the gap between the Hayashi track and the 
observed birth line of T Ta uri st ars in the Hertzsprung- 
Russel diagram (e.g Stahler 1988), it may not be as justi- 
fied in the infrared (IR). Including these embedded stars 
would however require substantial effort, since their SEDs 
are broa der than a single temperature blackbody (Shu et 
al. 1987 ) and may not be well-represented by our current 
library of stellar atmospheres. 

The correction for the parental cocoon only applies 
to the pre-main sequence tracks used in the range 0.8 < 
M/Mq < 7, and not to the low-mass tracks 0.08 < 
M/Mq < 0.8 (covering the pre-main sequence as well as 
later stages) due to the lack of detailed pre-main sequence 
time-scales for these objects. 



2.3. Horizontal branch morphologies 

Since the process responsible for mass-loss on the red gi- 
ant branch (RGB) is not properly understood, HB models 
are still highly uncertain. This becomes especially appar- 
ent in low metallicity globular clusters (GCs), where stan- 
dard models are unable to reproduce the wide range of HB 
morphologies observed. In order to set the stage for our 
approach to this problem, we present a short review of the 
on-going debate about the second parameter (metallicity 
being the first) responsible for determining the morphol- 
ogy- 

By assuming that more massive stars experience more 
intense mass loss along the RGB, age becomes a natu- 
ral cand idate for the second parameter (e.g. Fusi-Pecci & 
Renzini 1975 ). However, as has been argued by Buonanno 
et al. ( 1985 ) several parameters may play more or less 



significant role both within a specific cluster and be- 
tween clusters. As was discussed by Catalan & de Freitas 
Pacheco (1995), the two clusters M3 and M13 which 
have a similar age (Catelan & de Freitas Pacheco |1995 



Johnson & Bolte [199$ VandenBerg [1999| ; Grundahl [l999| j 



still show a large difference in HB morphology. Similarly, 
Catelan Q2000| ) and VandenBerg Q2000P , found that the 
difference in HB morphology between other galactic GCs 
proved too large to be accounted for by age dependent 
mass loss alone and that age could not be the dominant 
secondary parameter. 

It has been known for a long time that the Sculptor 
dwarf galaxy contains an unusually red HB p opulation 
for it s high a ge an d low metallicity (Zinn 1980 DaCosta 
In their recent study of the galaxy 



1984; Mateo 



1998) 



Hurley-Keller et al. ( 1999] ) note that while there is a strong 



concentration of these red HB stars towards the centre, no 
gradient in age is found. 

There seems to be a systematic tendency for GCs 
with blue extensions of the HB to contain super oxygen- 
poor stars, as measured by the [ O/Fe ] relative abundances 
(Catelan & de Freitas Pacheco [2000 ; Carretta & Gratton 
|1996| ). Since stars with these chemical properties appear 
to be losing mass at an abnormally high rate, affecting 
RGB evolution as well as later stages, detailed differences 
in chemical composition may well contribute to the second 
parameter effect. It is therefore necessary to understand 
the mechanism that influences the variations of chemi- 
cal composition of the stellar atmospheres. Cottrell & da 
Costa ( |1981| ) suggested a primordial origin of the varia- 
tions, possibly originating in the processed envelopes of 
intermediate-mass AGB stars that polluted the nearby 
cluster environment by different amounts. In this case the 
chemical variations of the cluster RGB stars would be of 



a purely stochastical nature. Peterson ( 1985 ) suggested 
that the second parameter might be related to stellar ro- 
tation. If so, such an effect may also explain the abundance 
anomalies through the induction of deep mixing processes. 
Most likely there would be a combination of primordial 
and induced origin. 

Another possibility is that the RGB stars are spun up 
by stellar companions resulting in an increased mass loss 
and a subsequent tendency to end up in the blue part of 
the HB. Rich (1997), however, regard the influence from 
binary companions to be insufficient. More importantly, 
and possibly correlated with the relative frequency of bi- 
naries and thus complementary, may be the influence of 
massive planets as argued by Soker & Harpaz (2000). 

The conclusion must be that there may be a rich 
spectrum of factors behind the second parameter effect 
and that age probably plays a minor role. Even so, most 
SEMs generally employ evolutionary tracks from either 
the Geneva or Padova group without further modifica- 
tions, and therefore implicitly assume the second parame- 
ter to be age. To study how this hypothesis may affect the 
predicted spectral evolution of galaxies we have chosen to 
adopt an opposing view. 
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In order to properly cover the whole range of observed 
morphologies at low metallicities, the theoretical HB-AGB 
tracks by Castellani et al. (1991) have been compared to 
the morphological classification of HBs for a large num- 
ber of Galactic GCs (Harris p9^ ). For Z=0.001, 0.004 
and 0.008, all tracks crossing the instability strip in a way 
that does not seem to be supported by these observations 
have been excluded from the set. Every star with initial 
mass below 1.7 M is then assumed to be equally likely 
to follow any of the remaining HB-AGB tracks. This pro- 
cedure will simulate a second parameter that is stochastic 
in nature and result in an age-independent distribution of 
stars along the HB that we claim will give a better cover- 
age of the many HB morphologies (especially at Z=0.001) 
confirmed observationally. 

Due to the scarce observational material on higher 
metallicity GCs, this technique breaks down at Z=0.020 
and Z=0.040, for which we have chosen to adopt the 
more mainstream approach provided by the use of Geneva 
tracks (Charbonnel et al. 199E). This should not result in 
any great inconsistency in terms of morphology, unless 
the HB becomes progressively bluer at metallicities above 
Z=0.008. 

Since our HB strategy assumes the existence of a large 
variety of horizontal branch morphologies within each sys- 
tem at low metallicities, caution must be taken when ap- 
plying the model to very small stellar systems (e.g. in- 
dividual globular clusters). It should of course be noted 
that if the second parameter truly is stochastic in nature, 
the uncertainties introduced by doing so may perhaps be 
no larger than those introduced when using conventional 
models. 

2.4. Theoretical stellar atmospheres 

The compilation of theoretical stellar atmospheres by 
Lejeune et al. ( [1998 ) has been used for stars with effec- 
tive temperature < 40000 K. These cover effective tem- 
peratures as low as 2000 K, surface gravities in the range 
-1.02 < log(g) < 5.5 and metallicities -3.5 < [A//H] < 
1.0. 

In order to match the five metallicities covered by our 
stellar evolutionary tracks, the original grid has been lin- 
early interpolated in \og(Z). For effective temperatures 
over 40000 K, th e non-LTE atmospheres by Clegg & 
Middlemass (1987]) have been employed after linear inter- 
polation between data points to yield the same resolution 
as for the Lejeune compilation (1221 data points between 
wavelengths 9.1-160000 nm). 

The two sets of theoretical atmospheres are not com- 
pletely compatible, and making the switch at 40000 K will 
cause an unavoidable discontinuity in the ultraviolet (UV) 
spectral evolution. This will manifest itself in a sudden 
drop of several orders of magnitude in flux at wavelengths 
below 23 nm as the stars of the stellar population cool off 
sufficiently to fall below this temperature. 



2.5. The nebular component 

The inclusion of a nebular component is of vital impor- 
tance for modelling regions of active star formation. Even 
so, many recent models choose to neglect the effect of ion- 
ized gas or to treat it onl y sup erficially. The Starburst99 
model of Leitherer et al. ( 1999 ) includes the nebular con- 
tinuum, but no emission lines. The PEGASE2 model of 



Fioc and Rocca-Volmerange ( 199E ) predicts a gas contin- 
uum as well as emission lines, but only by resorting to 
a nebular component pre-calculated for a set of fixed gas 
parameters 1 . 

Our nebular component (lines and continuum) is mod- 
elled using the predicted SED of the stellar population as 
a heat source for the photoionization code Cloudy ver- 



sion 90.05 (Ferland 199£). This makes our gas prescrip- 
tion far more flexible than existing models, and allows the 
possible impact of variations in the physical conditions of 
the gas on spectral evolution to be evaluated. Variable in- 
put parameters include the covering factor, filling factor, 
gas metallicity and hydrogen number density. The cover- 
ing factor simulates the existence of possible holes in the 
nebula, whereas the filling factor quantifies the degree of 
graininess of the gas. 

When applying Cloudy, some important approxima- 
tions have been used: The stars are assumed not to be 
spatially mixed with the gas, the hydrogen density is as- 
sumed to be constant in time, the nebula to be ionization 
bounded and the geometry to be spherical. 

2.6. Numerical method 



Following Fioc & Rocca-Volmerange ( 1997 ), the 
monochromatic flux of a galaxy at time i ga i may be 
written: 

/■tgal /"». 

F\(t g3 x,T, a) = / / SFR(t ga i - i star , t)$(to, a) ■ 

■f\(m,t stllr ) dm di s tar 



Jm 



provided that $(m, a) is normalized to 1 Mq. Here, 
fxi^m^tst^) represents the monochromatic flux of a star 
with initial mass m at wavelength A and age i s t a r- 

In order to perform spectral evolutionary synthesis, 
at least some part of this expression is commonly dis- 
cretized. Discrctizing both integrals leads to a somewhat 
outdated method of evolutionary synthesis, whose draw- 
backs were thoroughly investigated by Chariot & Bruzual 
( 1991 ) , whereas discretization of only one integral leads to 
the more modern isochrone (time integral discretized) and 
isomass (mass integral discretized) methods. This model 
is based on the latter technique, where the expression for 
F\ (t ga i, t, a) is rewritten as: 

p-i 
F\(t gah T, a) =y^ SFR(i gal - t sta r,i, r) ■ 

i=\ 



1 However, a coupling of the two codes PEGA SE2 and 
Cloudy was recently published by Moy et al. ( 2001 ) and ap- 
plied to starbursts and HII galaxies. 
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^$(mj,a)(m j+ i -my) 



3 = 1 



f\(mj,t sta ,r) di s 



Here £ sta r,i = 0, £ s tar, P = igai, mi = mi and m q = m u . The 
number of time and mass intervals used in the evaluation 
is represented by p — 1 and q — 1, respectively. 

The problem of performing a simultaneous discretiza- 
tion of both time and mass integrals is that it may cause 
oscillating colours and luminosities when modelling short 
bursts of star formation. These phenomena arise when 
stars on one track abruptly move into a different evolu- 
tionary phase that stars on the consecutive track will not 
reach until later. Oscillatio ns of t his kind were a major rea- 
son for Chariot & Bruzual ( 1991 ) to propose the technique 
of isochrone synthesis as a better alternative, although the 
isomass method should theoretically be just as valid. This 
was pointed out by Fioc & Rocca-Volmerange ( 1997 ) who 
also claimed to have tested both methods with identical 
results. There is, however, always the risk that hidden ap- 
proximations present in the implementation of these algo- 
rithms may have an impact under certain circumstances. 
Since our ambition has been to develop an independent 
code that allows critical evaluation of the uncertainties 
involved in the predictions, we have therefore preferred 
the use of the isomass method to the more mainstream 
isochrone approach. 

In this technique, it is necessary to ensure that mj+i- 
rrij is sufficiently small so that equivalent evolutionary 
phases of consecutive masses overlap. For this reason, the 
original mass grid is interpolated to yield a set of stellar 
evolutionary tracks appropriate for each metallicity and 
set of time steps. 

3. The impact of nebular emission 

We here investigate the magnitude and duration of nebu- 
lar effects on broadband colours, and test the sensitivity of 
colours and emission line equivalent widths to the metal- 
licity and physical conditions of the gas. 

In Table 1, the duration of observable nebular effects in 
different photometric filters is calculated for two different 
metallicities and three star formation scenarios. As seen, 
inclusion of the nebular component may alter the photo- 
metric properties of galaxies from the UV to the near-IR 
during substantial periods of time. These time-scales may 
serve as estimates to the age of a galaxy after which neb- 
ular effects may be neglected, a point often overlooked in 
the analysis of high-redshift objects. For each filter, the 
time-scale of significant nebular emissions is determined 
by the temporal evolution of the bright emission lines and 
the relative importance of continuum emission. For in- 
stance, the sharp decline in the time-scales between filters 
R and I in the most prolonged star formation scenario is 
caused by the presence of Ha in the R-band and the lack 
of durable lines of sufficient luminosity in I. The subse- 
quent increase in time-scales towards longer wavelengths 
is due to increasing continuum emission throughout filters 
I, J and K. 



Table 1. Duration of measurable contributions in Johnson 
filters from nebular emission. The contributions have been 
evaluated at metallicities Z=0.001 and Z=0.020 for three 
different star formation scenarios: A short burst of con- 
stant star formation (r=100 Myr) (S) and exponentially 
decreasing SFRs with e-folding decay rates r=l Gyr (El) 
and r=14 Gyr (E14). All scenarios were computed us- 
ing a Salpeter IMF over the mass range 0.08-120 Mq, 
10 10 Mq available for star formation, Z s t als = Z gas , cover- 
ing factor 1.0, filling factor 1.0 and n(H) = 100 cm -3 . The 
entries r (filter) measure the time in Myr after which the 
discrepancy between models with and without the nebular 
component become negligible, i.e. smaller than the typical 
observational uncertainty of 0.05 magnitudes. 
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By varying the parameters regulating the physical 
properties of the gas we have also investigated the sen- 
sitivity of the predicted colours to such factors. The sen- 
sitivity to reasonable variations of the covering factor, fill- 
ing factor and hydrogen number density for a short burst 
scenario is illustrated in Fig. EL As expected, the most dra- 
matic effects come from variations in the covering factor, 
shifting the colours towards those of the a naked stellar 
component for low parameter values. Although the effec- 
tive value of this parameter remains unknown, it is fair to 
say that no studies to date have unraveled large Lyman 
continuum escape fractions (which would be one possi- 
ble signature of significant departures from covering factor 
unity) from galaxies in the local univcrsen. More trouble- 
some is perhaps the fact that variations in filling factor 
and hydrogen density may alter the colours by as much 
as 0.1 magnitudes. The dependence of predicted colours 
on gas metallicity, when the stellar metallicity is held con- 
stant, is plotted in Fig. g. Here, the dependence is shown 
to be even stronger than for the covering factor, peaking at 
0.8 magnitudes in V — R. Since even modest variations of 
the physical properties of the gas may have such a signifi- 
cant impact on the predictions, the common procedure of 
assuming parameter values typical of star-forming regions 
without a more careful analysis of the possible parameter 
space should be considered quite a hazardous exercise. 

In Fig. y[ we present the evolution of the Ha, H/3 and 
[OIII]A5007 line equivalent widths for different combina- 
tions of gas parameters and metallicities. When the stel- 



2 But see the recent study by Steidel et al. ( 2001 ) for a dif- 
ferent picture at high redshift. 
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Fig. 1. The evolution of the U - B, B - V, V - R, V - I, 
V — J and V — K colours predicted by our model at 
Z=0.001 when the covering factor, filling factor and hy- 
drogen number density (cm -3 ) are allowed to take on 
the following values: 1.0, 1.0, 100 (solid); 0.5, 1.0, 100 
(dashed); 1.0, 0.1, 100 (dash-dotted); 1.0, 1.0, 10 (dotted); 
0.5, 0.1, 10 (solid with dots). All evolutionary sequences 
assume a Salpeter IMF (0.08-120 M ), a short burst of 
constant star formation (t=100 Myr), 10 10 M Q available 
for star formation and Z s t a , rs =Z f r as . 



lar and nebular metallicities are considered to have the 
same value, the Ha and H/3 equivalent widths only show 
a modest metallicity dependence, whereas the metallic- 
ity dependence of [OIIIJA5007 is rather pronounced. In 
the case when the stellar metallicity is held constant at 
Zstars=0.001 and the gas metallicity Z gas is allowed to 
vary, both Ha and H/3 equivalent widths increase with in- 
creasing Z gas . For the [OIII]A5007 line, the highest equiva- 
lent widths are instead produced by the intermediate Z gas . 
The Ha and H/3 equivalent widths are apparently only 
sensitive to the covering factor, whereas the [OIII]A5007 
line also shows a slight dependence on filling factor and 
density. 

Stasihska & Leitherer ( 1996] ) investigated the evolution 
of the [OIIIJA5007 line equivalent width in a similar way, 
predicting monotonically decreasing equivalent widths for 
lower metallicities, but oscillating and increasing widths 
at solar metallicity, thereby hampering the use of this line 
as an age indicator for high-metallicity objects. This be- 
havior is not reproduced by our model, where the width 
shows a monotonic decrease with age for all metallicities, 
although not a very strong one at Z=0.040. The origin 
of this discrepancy is not known, but could possibly be 
connected to the lack of realistic atmosphere spectra for 
Wolf-Rayet stars in the library of Clegg & Middlemass 
( P87p - 




Fig. 2. The evolution of U -B, B- V, V- R, V-I, V - J 

and V — K colours predicted by our model for a stellar 
metallicity Z stars =0.001 and gas metallicities Z gas =0.001 
(solid), Z gas =0.004 (dashed), Z gas =0.008 (dash-dotted), 
Z gas =0.020 (dotted), Z gas =0.040 (solid with dots), assum- 
ing covering factor 1.0, filling factor 1.0, hydrogen number 
density 100 cm~ 3 and the same IMF and star formation 
history as in Fig. H. 



4. The impact of pre-main sequence evolution 



Despite sporadic evidence (e.g. Chariot 1996) pointing in 
a different direction, the PMS is normally assumed to be 
too short-lived or faint to seriously affect the integrated 
SED of stellar populations. By comparing evolutionary 
scenarios with and without inclusion of this phase, we have 
tested the validity of this assumption. The results are pre- 
sented in Tables 2 and 3, where the contribution to the 
integ rated luminosity from the PMS tracks by Bernasconi 
1996J , subject to the approximations specified in section 
2.2, has been evaluated in different photometric filters. 
Table 2 compares scenarios which invoke these tracks to 
cover the time gap in the Geneva set from the Hayashi- 
track to the ZAMS, with and without the inclusion of the 
light emitted during this phase. Table 3 instead compares 
scenarios which include the light emitted and time spent 
on the Bernasconi tracks to those where all stars are as- 
sumed to climb the ZAMS at time equal to zero. Which 
of these tables that best describe the PMS contribution 
to the SED depends on how the time zero is defined in its 
absence. 

Both tables evaluate the PMS contributions at metal- 
licities Z=0.001 and Z=0.020 for two different IMFs: A 
standard scenario computed using a Salpeter IMF over 
the mass range 0.08-120 Mq, and an extreme scenario, de- 
signed to maximize the effects of the lowest mass stars us- 
ing a different IMF (a=-2.85) and the smaller mass range 
0.08-25 Mq. Both scenarios were calculated using a short 
burst of star formation (r=100 Myr) and time steps of 1 
Myr. 
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Table 4. Time-weighted averages of the colour discrepancies between the evolutionary sequences from our model 
(Z2001), PEGASE2 (PEG2) and BC96 seen in the comparison of Fig. H. Each column represents the comparison 
between two sequences. 
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Fig. 3. The evolution of the Ha, H/3 and [OIII]A5007 
equivalent widths for the same IMF and star forma- 
tion history as in Fig. [y. Left column: Metallicities 
(Z stars =Z gas ) Z=0.001 (solid), Z=0.004 (dashed), 
Z=0.008 (dash-dotted), Z=0.020 (dotted), Z=0.040 
(solid with dots). Middle column: Stellar metal- 
licity Z st ars=0-001 and gas metallicities Z gas =0.001 
(solid), Z gas =0.004 (dashed), Z gas =0.008 (dash-dotted), 
Z gas =0.020 (dotted), Z gas =0.040 (solid with dots), Right 
column: Metallicity (Z stllIS =Z gas ) Z=0.001, covering 
factor, filling factor and hydrogen number density (cm ) 
allowed to take on the following values: 1.0, 1.0, 100 
(solid); 0.5, 1.0, 100 (dashed); 1.0, 0.1, 100 (dash-dotted); 
1.0, 1.0, 10 (dotted); 0.5, 0.1, 10 (solid with dots). 



Regardless of evaluation method, the PMS can be seen 
to produce measurable signatures in the near-IR evolu- 
tion during the first few Myr of a standard IMF, short 
burst scenario. If more extreme IMFs are to be modelled, 
the effects become significant in all filters and for much 
more prolonged periods of time. The PMS should there- 
fore certainly not be neglected when very young stellar 
populations are being modelled. 

5. Comparison to other models 

To illustrate the uncertainties still present in contempo- 
rary spectral evolutionary synthesis, we present a compar- 



Table 2. Duration and amplitude of measurable contribu- 
tions in Johnson filters from the PMS tracks by Bernasconi 



1996, evaluated by comparing evolutionary sequences with 



and without the inclusion of the light emitted during this 
stage. The entries max(A(filter)) measure the maximum 
difference in magnitudes between models with and with- 
out PMS light (negative max(A(filter)) indicate an in- 
crease in luminosity due to the inclusion of light emitted 
during this stage). The entries r (filter) measure the time 
in Myr after which the discrepancy between models with 
and without pre-main sequence light become negligible, 
e.g. smaller than the typical observational uncertainty of 
0.05 magnitudes. 
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ison between colour predictions from our model and those 
produced by the widely-used codes of Bruzual & Chariot 



( 1996 ) (hereafter BC96) and Fioc & Rocca-Volmerange 
( 1999 ) (PEGASE2) when similar input parameters are be- 
ing used. Johnson filters have been employed throughout. 

PEGASE2 enables a comparison between stellar and 
nebular components at all five metallicities used by our 
model, whereas BC96 only allows comparison between 
stellar components at Z=0.004, 0.008 and 0.020. 

In Fig. the photometric evolution predicted at 
Z=0.004, 0.008 and 0.008 for a pure stellar population 
undergoing a short burst of star formation is plotted for 
all three codes. This scenario is more suited for a quan- 
titative comparison than any of the more continuous star 
formation histories, since a larger range of masses will be 
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Table 3. Duration and amplitude of measurable contribu- 
tions in Johnson filters from the PMS tracks by Bernasconi 



1996, evaluated by comparing evolutionary sequences in- 



cluding these tracks to those assuming all stars to start on 
the ZAMS t—0. The entries max( A (filter)) measure the 
maximum difference in magnitudes between models with 
and without PMS tracks (negative max( A (filter)) indi- 
cate an increase in luminosity due to the inclusion of this 
phase). Otherwise same as Table 2. 
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able to contribute to the integrated colours. Fig. |5| con- 
tains the same comparison at 2=0.001 and Z=0.QAQ for 
our code and PEGASE2. 

Our code, PEGASE and BC96 all seem to agree fairly 
well in U — B and B — V colours, with differences peaking 
at 0.2 magnitudes. However, in V — K the differences may 
reach as high as 0.7 magnitudes. Table 4 lists the time- 
averaged discrepancies between the three models along 
the colour sequences of Fig. [|. BC96 and PEGASE2 gen- 
erally seem to be in better agreement with each other 
than with our model. This is hardly surprising, since our 
model is mainly based on tracks from the Geneva set, 
whereas BC96 and PEGASE2 both employ stellar evo- 
lutionary tracks from the Padova set. Only in U — B does 
the mean difference between our model and PEGASE2 
become smaller than that of BC96 and PEGASE2. Apart 
from the Geneva-Padova discrepancy, the pronounced dif- 
ferences in V — K between our model and the other two 
are enhanced by our use of evolutionary stages and strate- 
gies not adopted in these codes. The inclusion of PMS 
makes V — K redder during the first few Myr of evo- 
lution, whereas the bluer V — K seen in the age inter- 
val 0.1-1 Gyr is likely due to the different prescriptions 



used for the TP-AGB (Fioc & Rocca-Volmerange |l997j). 
Finally, the effects of our extended HB morphologies at the 
three lowest metallicities become important after about 1 
Gyr, producing significantly bluer V — K at Z=0.004 and 
Z=0.008. The slightly redder V-K predicted at Z=0.001 
indicates a decreasing sensitivity to the HB morphology at 
this metallicity and the dominance of other evolutionary 



stages. Unfortunately, the discrepancies stemming from 
the different approaches to the TP-AGB and second pa- 
rameter problem can not easily be disentangled from the 
different stellar tracks, atmospheres and computational 
methods used. To properly resolve this issue would re- 
quire a much more careful analysis that is not within the 
scope of this paper. 

In Fig. a we compare the photometric evolution pre- 
dicted for stellar populations experiencing a short, con- 
stant burst of star formation at Z=0.001 and Z=0.020 
using our code and PEGASE2, with and without the in- 
clusion of the nebular component. Since the input param- 
eters used by our photoionization model cannot easily be 
translated into the language of PEGASE2, nebular contri- 
bution in PEGASE2 was computed using default settings, 
while we preferred the use of n(H)=100 cm -3 , a filling fac- 
tor of 1.0, a covering factor of unity and 10 10 Mq available 
for star formation. Even though both models start out 
with identical metallicities (Z ststIS = Z gas ) , evolution in the 
gaseous metallicity of PEGASE2 cannot not be prevented, 
making the comparison less adequate at high ages. In the 
evolutionary sequences used here, PEGASE2 will reach 
Z gas =0.004 at t w 95 Myr in the Z stars =0.001 scenario 
and Z gas =0.040 at t « 75 Myr when Z sta rs=0.020. 

As seen in Fig. pi properly taking the effects of nebu- 
lar emission into account becomes crucial during the star- 
forming phase. The differences in colour between scenar- 
ios with and without nebular emission peak at 0.8 (0.2) 
in 17 - £ for our code (PEGASE2), 0.7 (0.6) in B - V 
and 0.7 (1.3) in V — K. Due to the different input pa- 
rameters, stellar SEDs and photoionization models used 
the final colours are however very different, especially in 
U — B, where the sequences actually shift in different di- 
rections when nebular emission is added. The significantly 
bluer U — B produced by PEGASE2 appears to stem from 
a much higher luminosity predicted for the [OIIJA3727 
emission line throughout the entire star formation episode 
([OII]A3727/H/3 « 3 already at t=\ Myr). Even though a 
direct comparison of these colours to those seen in star- 
burst galaxies may be complicated by the possible pres- 
ence of dust, it should be noted that colours as blue as the 
ones predicted by PEGASE2 are not typically observed. 
Out of the ~ 30 blue compact galaxies (assumed t o have 
U-B < -0.1) in the sample of Bergvall & Olofsson fll986| ) 
and 17 blue compact dwarf galaxies (for which this colour 
was available) in Thuan (1983), none have U — B < — 1. 

Since reliable tests of SEMs are yet to be developed, 
there is no simple way around these model discrepan- 
cies. Our recommendation would therefore be to conduct 
multi-model consistency checks of all important conclu- 
sions whenever comparable codes are available. 

6. Comparison to observations 

To check the performance of our model in confrontation 
with reality, we have made a comparison of predicted 
broadband colours to those of observed globular clus- 
ters (GCs) and low-surface brightness galaxies (LSBG). 
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Fig. 4. The evolution of the U - B, B - V and V - K 
colours predicted by our code (solid), BC96 (dashed) and 
PEGASE2 (dotted) at Z=0.004, 0.008 and 0.020. All mod- 
els have been computed assuming a Salpeter IMF (0.08- 
120 M Q for our code, 0.1-120 M© for PEGASE2, 0.1-125 
Mq for BC96), a short burst of constant star formation 
(t=100 Myr) and no nebular contribution. 



Predicted line ratios and equivalent widths have been com- 
pared to observed HII galaxies. Unless otherwise stated, 
all models assume Z &as =Z sia , TS , covering factor 1.0, filling 




Fig. 5. The evolution of the U - B, B - V and V - K 

colours predicted by our code (solid) and PEGASE2 (dot- 
ted) at Z=0.001 and 0.040. Otherwise same as Fig. §. 



and a Salpeter IMF (0.08-120 



factor 1.0, n(H)=100 cm" 

M Q ). 

6.1. Broadband colours 



Even though our model, due to the statistical treatment 
of HB morphologies, cannot be expected to reproduce 
the average colours of low-metallicity globular clusters, 
as discussed in section 2.3, we would expect the predicted 
colours to fall somewhere within the range observed. In 
order to test this property, we have compared our evo- 
lutionary sequences to the large set of GCs in M31, as 
compiled by Barmby et al. ( |2000| ). Since the metallicities 
and dereddened colours derived for theses objects are only 
available in figures, the comparison is carried out with re- 
spect to the approximate outer boundary of the region 
inhabited by the GCs the B — V, V — K two-colour dia- 
gram, especially suited to test the treatment of late stages 
of stellar evolution. The boundaries have been subjected 
to metallicity binning in V — K but not in B — V . Since 
our main interest is to make sure that the extended HB 
morphologies produced by our approach to the second pa- 
rameter problem do not produce sequences that are too 



10 



Zackrisson et al.: Spectral galaxy evolution including the effects of nebular emission 




Fig. 6. The evolution of the U - B, B - V, V - K 
colours predicted with/ without the inclusion of the nebu- 
lar component by our code (solid/dashed) and PEGASE2 
(dotted/dash-dotted) at Z=0.001 and Z=0.020, assum- 
ing a short burst of constant SFR (r=100 Myr), a 
Salpeter IMF (0.08-120 M Q for our code, 0.1-120 M for 
PEGASE). 



blue in V — K at high ages, this is of no major concern. 
As seen in Fig. J7|, the evolutionary sequences at Z=0.001, 
0.004 and 0.008 all fall well within the range of colours 
seen in this sample for ages higher than about 5 Gyr. 
We now turn our attention to a set of blue LSBGs 



(Bergvall et al. |1999| ; Ronnback & Bergvall |1994J) which 
despite being actively star-forming and metal-poor are 
believed to have ages higher than 2 Gyr (Bergvall & 
Ronnback 1994). Standard models of chemical evolution 
cannot easily account for this phenomenon without fine- 
tuning the unknown rates of infall to and mass-loss from 
these galaxies. A constant metallicity model may there- 
fore be just as appropriate for describing their evolution. 
A detailed analysis of the star formation history of these 
objects is however not within the scope of this paper. 
Here, we are simply want to establish that the colours ob- 
served can be reasonably well reproduced by our model. 
The low level of dust extinction observed in these objects 
allows the predicted colours to be directly compared to 




Fig. 7. Regions in B — V vs V — K space inhabited 
by confirmed GCs in M31, subject to metallicity bin- 
ning in V — K . Left: Evolutionary sequence (solid) pre- 
dicted at Z=0.001 and the region inhabited by objects 
in the metallicity range Z=0. 0005-0. 002 (shaded area). 
Middle: Evolutionary sequence at Z— 0.004 and objects in 
the range Z=0. 002-0. 006. Right: Evolutionary sequence 
at Z=0.008 and objects in the range Z=0.006-0.014. All 
models assume a Salpeter IMF (0.08-120 Mq, a short 
burst of constant star formation (r=10 Myr) and 10 5 M Q 
available for star formation. For each model sequence, cir- 
cles mark the predicted colours at 1 Gyr, diamonds at 5 
Gyr, squares at 10 Gyr and pentagrams at 15 Gyr. 

the observations without any correction for reddening. In 
Fig. [§] we compare the location of the observed LSBGs 
in the B — V, V — J two-colour diagram to the evolu- 
tionary sequences predicted for scenarios believed to span 
the likely range of present star formation in this sample. 
With the exception of one object, separated from the oth- 
ers in V — J, the model appears quite compatible with the 
colours observed. Even though the scenario with shorter 
time-scale of star formation probably underestimates the 
current SFR in these objects, it places firm constraints 
on the lowest ages allowed, indicating - in agreement with 
previous findings - that these galaxies cannot be young. 

6.2. Line ratios and equivalent widths 

For the spectroscopic data, we have been usi ng th e com- 
pilation of HII galaxies by Masegosa et al. ( 1994 ), con- 
sisting of 121 objects, all corrected for intrinsic reddening. 
This data-set was also analyzed by Stasinska & Leitherer 
( 1996J ) , which allows the performance of their model to be 
directly compared to ours. 

Since the starburst regions sampled in this data-set 
are of varying spatial extent, but typically smaller than 
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Fig. 8. B - V vs V - J for 9 blue LSBGs (dots). Lines 
are evolutionary sequences predicted at .£=0.001 for expo- 
nentially declining SFRs with r=l Gyr (solid) and r=14 
Gyr (dashed), assuming 10 10 M Q available for star forma- 
tion. Along each sequence, ages of 1, 5, 10 and 15 Gyr have 
been illustrated with markers (circles for solid, squares for 
dashed). 

the host HII galaxy, the total mass available for star for- 
mation has been considered a free parameter. We have 
furthermore assumed that the contribution to the SED 
from the older, underlying stellar component is insignifi- 
cant. 

Fig. |9j presents observations and evolutionary se- 
quences in two diagnostic diagrams, [OIIIJA5007/H/3 vs 
W(H/3) and [OII]A3727/H/3 vs W{R/3). The evolutionary 
sequences of [OHI]A5007/H/3 vs W(H/3) predicted for star- 
burst masses in the range 10 5 — 10 10 Mq all seem to 
agree reasonably well with the observations. The value 
of W(T3.P) at which the predicted [OIII]A5007/H/? starts 
to drop is sensitive to the duration of the starburst, and 
the value of [OIILJA5007/H/3 for objects located at VF(H/?) 
lower than 5 nm may indicate more prolonged star forma- 
tion scenarios than the one assumed here. 

The [OII]A3727/H/3 vs W(R(3) diagram presents a 
much more troublesome picture. Neither the slope nor the 
absolute values of the observed [OII]A3727/H/9-sequence 
can easily be reproduced by our model. The line ratio is 
significantly underpredicted during the entire duration of 
the starburst for all but the lowest mass model, which 
performs poorly in [OIIIJA5007/H/3 vs W(H/3), especially 
in the lower metallicity bin. A similar trend is seen in 
Stasinska & Leitherer, where only the 10 3 Mq model - also 
rating poorly in [OIII]A5007/H/3 vs W(R/3) - was able to 
reproduce the high values of [OII]A3727/H/3. Assuming 
masses and efficiencies typical of star-forming regions, 
these low-mass models appear quite unlikely. It seems rea- 
sonable to believe that the stellar mass in this sample 
should be at least 10 5 -10 6 Mq. In Fig. [l(| we therefore in- 
vestigate the dependence of the 10 5 M Q model on different 
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Fig. 9. Evolutionary sequences of [OIIIJA5007/H/3 vs 
W{H/3) and [OII]A3727/H/3 vs W(H/3) assuming 10 3 M© 
(solid), 10 5 Mq (dashed), 10 8 M Q (dash-dotted) and 
10 10 Mq (dotted) available for star formation. All models 
use a Salpeter IMF (0.08-120 Mq) and a short burst of 
constant SFR (r=10 Myr). The observed sample of HII 
galaxies is represented by dots. Asterisks mark ages of 
0.5, 5, 10 and 15 Myr. Left: Evolutionary sequences as- 
suming Z = 0.001 and observed metallicities in the range 
Z=0. 0005-0. 002. Right: Evolutionary sequences assuming 
2=0.004 and observed metallicities in the range Z=0.002- 
0.006. 



gas parameters, and show that the discrepancies between 
the performance in the two diagnostic diagrams cannot 
easily be resolved this way. Since this curiosity has now 
been confirmed by two independent codes, using different 
stellar evolutionary tracks, stellar atmospheres and pho- 
toionization models, this may indicate something missing 
from our understanding of the physics of these objects. 
Other observat ional studies of emission line galaxies (e.g. 
Kniazev et al. 2001 ) show similar trends in these diag- 
nostic diagrams, indicating that severe problems with the 
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actual data-set can probably be ruled out. A plausible ex- 
planation for the discrepancy could however be the exis- 
tence of spatial variations in the physical conditions (e.g. 
density and filling factor) between the [OIII] and [Oil] 
regions of the nebula. There could also be a mixture of 
several distinct HII regions with different characteristics 
sampled by the slit of the spectrograph, thereby making 
the model assumption of a single, homogeneous HII re- 
gion inappropriate. Part of the nebular emission may fur- 
thermore originate from regions ionized by shocks or non- 
thermal sources. As a test of this last possibility, we have 
compared the observed emission line ratios to those pre- 
dicted by a simple model using a mixture of photoionized, 
shocked and power-law ionized regions. Three regions were 
mixed, either 3 HII, 2 HII + 1 shock, 2 HII+1 power-law 
or 1 HII + 1 shock + 1 power-law region. Det ails ab out 
the procedure may be found in Bergvall et al. ( 1986 ). In 
« 20% of objects, combinations including power-laws or 
shocks gave significantly better fits than did pure mixtures 
of HII regions. This indicates that shocks and non-thermal 
sources could, to some degree, influence line ratios in this 
sample and explain part of the consistency problems be- 
tween models and observations. 

7. Summary and conclusions 

We have presented the first version of a new spectral evo- 
lutionary code for galaxies in the age range 0-15 Gyr. The 
model represents an improvement with regard to existing 
codes in the inclusion of PMS evolution and a sophisti- 
cated nebular component. 

The time durations of significant nebular emission 
in different photometric filters are estimated for several 
star formation scenarios and metallicities. These estimates 
may serve as a guide to the age at which the common 
assumption of negligible nebular emission in early-type 
galaxies breaks down. The sensitivity of the colours pre- 
dicted to the physical properties of gas (i.e hydrogen den- 
sity, metallicity, covering and filling factors) is also tested, 
revealing a high sensitivity to the metallicity and covering 
factor, and a lower, yet observationally significant, sensi- 
tivity to the density and filling factor assumed. This has 
implications for the modelling and analysis of young galax- 
ies, where normally gas parameter values only typical of 
star-forming regions are applied, and stresses the need for 
a more careful analysis of the possible parameter space. 

The predicted evolution of the Ha, H/3 and [OIII]A5007 
line equivalent widths is presented for different metallici- 
ties and values of gas parameters. The Ha, H/3 equivalent 
widths are found only to be sensitive to the covering fac- 
tor whereas [OIIIJA5007 appears sensitive to the choice 
of other parameters as well. It is noted that the oscil- 
lating and increasing equivalent width of [OIII]A5007 pre- 
dicted by Stasifiska & Leitherer ( 1996 ) at solar metallicity, 
thereby hampering the use of this line as an age indicator, 
is not reproduced by our model. 

The contribution from the PMS phase to the inte- 
grated properties of a young stellar populations is eval- 
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Fig. 10. The gas parameter dependence of the 10 5 Mq 
model from Fig. g. The covering factor, filling factor and 
hydrogen number density (cm -3 ) were allowed to take on 
the following values: 0.5, 1.0, 100 (solid); 1.0, 0.1, 100 
(dashed); 1.0, 1.0, 10 (dash-dotted); 0.5, 0.1, 10 (dotted). 



uated and found to be substantial in the near-IR during 
the first few million years of a standard IMF starburst 
lasting 100 Myr, e.g. increasing the luminosity in the K- 
band by as much as 0.34 magnitudes at Z=0.001 and 0.16 
at Z=0.020 when compared to the case where all stars 
are assumed to simultaneously start on the ZAMS at time 
equal to zero. If more extreme IMFs are to be modelled, 
the effects become significant in all filters and for much 
more prolonged periods of time. 

A comparison to the colour predictions of the stel- 
lar components from two similar models, BC96 and 
PEGASE2, reveal a number of discrepancies between the 
codes, especially in the near-IR, where predictions are 
highly affected by uncertainties in the treatment of late 
stages of stellar evolution. Since reliable tests of SEMs 
are yet to be developed, this result emphasizes the need 
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for multi-model consistency checks whenever comparable 
codes are available. 

Finally, the colours predicted by our model are tested 
against observed samples of GCs and LSBGs, showing an 
acceptable agreement. Line ratios and equivalent widths 
are then compared to a set of HII galaxies. Even though 
the model performs well in the [OIII]A5007/H/3 vs W(B./3) 
diagnostic diagram, the difficulty in simultaneously fitting 
[OIII]A5007/H/3 vs W(Rj3) and [OII]A3727/H/3 vs W(H/3) 
noted in Stasiriska & Leitherer ( 1996| ) does however persist 
in our study as well. This may indicate something missing 
from our understanding of the physics of these objects. 

We conclude that our model is appropriate for the 
modelling of young as well as old galaxies in the low- to in- 
termediate mass range. Due to the large uncertainties still 
prevailing the field of spectral evolutionary synthesis (e.g. 
the near-IR evolution of stellar populations, the colours 
produced by inclusion of the nebular component), we do 
not recommend the use of any model by itself in the anal- 
ysis of galaxies, and stress that all important conclusions 
should be cross-checked with other codes as well. 
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